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We study the singly-Cabibbo-suppressed decays Aj — pn and AZ — pm? using AZ A; pairs 


produced by e*e^ 


collisions at a center-of-mass energy of \/s = 4.6 GeV. The data sample was 


collected by the BESIII detector at the BEPCII collider and corresponds to an integrated luminosity 
of 567 pb^!. We find the first evidence for the decay At — pn with a statistical significance of 4.20 
and measure its branching fraction to be B(A — pn) = (1.24 + 0.28(stat.) + 0.10(syst.)) x 107°. 
No significant AZ — pz? signal is observed. We set an upper limit on its branching fraction 
B(At — pn?) < 2.7 x 1074 at the 9096 confidence level. 


PACS numbers: 14.20.Lq, 13.30.Eg, 12.38.Qk 


Weak decays of charmed baryons provide a unique 
testing ground for different theoretical models and ap- 
proaches, e.g. the quark model approach to non-leptonic 
charm decays and Heavy Quark Effective Theory EHA. 
The charmed baryon ground state A? was first observed 
in 1979 [al [9], but, compared to the rapid advances of 
charmed mesons, progress in the studies of the charmed 
baryons has been relatively slow due to a lack of experi- 
mental data and the additional difficulties of three con- 
stituent quarks in theoretical calculation. The accuracy 
of At branching fractions (BFs) has long been poor for 
the Cabibbo-favoured (CF) decays, and even worse, with 
uncertainties at the 4096 level, for the singly-Cabibbo- 
suppressed (SCS) decays [19]. As a consequence, it is 
neither possible to test the BFs predicted by different 
theoretical models, nor to determine the effects of final- 
state interactions (FSI). It is therefore essential to im- 
prove the accuracy of these BFs for AZ decays and to 
search for new decay modes. The absolute BFs of twelve 
At CF hadronic decay modes have been measured by the 
BESIII collaboration with much improved precision ; 

The SCS decays At — pn and pz? have not yet 
been studied experimentally. These two decays pro- 
ceed predominantly through internal W-emission and 
W-exchange diagrams, which are non-factorizable and 


not subject to color and helicity suppression in charmed 
baryon decay. Some theoretical models [3 H [13] [13], 
predict the BFs of these two process under different as- 
sumptions (the flavor SU(3) symmetry, FSI) obtaining 
different results. Therefore, measurements of these BFs 
will help us to understand the underlying dynamics of 
charmed baryon decays and distinguish between the dif- 
ferent models. Furthermore, the ratio of BFs of these 
two decays, which is expected to be relatively insensitive 
to the values of input parameters in the theoretical cal- 
culation, is an excellent probe to distinguish between the 
different models. 

In this Letter, we present the first experimental in- 
vestigations of the SCS decays At — pn and pr°. We 
use a data sample of ete” collisions at a center-of-mass 
(c.m.) energy of ys = 4.6 GeV with an integrated 
luminosity of 567 pb-! collected by the BESIII 
detector at the BEPCII collider. Taking advantage of 
the excellent BESIII detector performance and the clean 
environment just above the mass threshold to produce 
ACA. pairs, a single-tag method, (i.e., reconstruction 
of only one A, in the At A; pairs) is used to increase 
the detection efficiency and acquire more A, candidates. 
'Throughout the text, the charge conjugate states are al- 
ways implied unless mentioned explicitly. 
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BESIII is a cylindrical spectrometer, consisting of 
a small-celled, Helium-based main drift chamber (MDC), 
a plastic scintillator Time-of-Flight system (TOF), a 
CsI(TI) electromagnetic calorimeter (EMC), a supercon- 
ducting solenoid providing a 1.0 T magnetic field, and a 
muon counter. The charged particle momentum resolu- 
tion is 0.596 at a transverse momentum of 1 GeV/c and 
the photon energy resolution in the EMC is 2.596 (596) 
in the barrel (endcap) region for 1 GeV photons. A more 
detailed description of the BESIII detector is given in 
Ref. [e]. 

High-statistics e*e^ annihilation Monte Carlo (MC) 
samples, generated by the GEANT4-based [18| MC 
simulation package BOOST [29], are used to investigate 
the backgrounds, to optimize the selection criteria, and 
to determine the detection efficiencies. The efe anni- 
hilation is simulated by the MC generator KKMC [21], 
taking into consideration the spread of the beam energy 
and the effect of the initial-state radiation (ISR). Inclu- 
sive MC samples, consisting of A* Az events, charmed 


meson D pair production, ISR returns to lower mass 


charmonium(-like) states, and continuum QED pro- 
cesses ete” — qq (q = u,d,s), are used to study the 
backgrounds. All known decay modes are generated with 
EVTGEN [23 [23] with BFs being the values of the Particle 
Data Group (PDG) [1d], and the remaining unknown de- 
cay modes are generated by LUNDCHARM . The signal 
MC samples of ete~ — ATA; are produced with one A, 
decaying to the final states of interest, py or pr°, and the 
other A, decaying generically to any of the possible final 
states. 

Charged tracks, reconstructed from hits in the MDC, 
are required to have a polar angle 0 satisfying | cos0| < 
0.93 and a point of closest approach to the interaction 
point within +10 cm along the beam direction (V?) and 
1 cm in the plane perpendicular to the beam (V,). In- 
formation from the TOF is combined with the ionization 
energy loss (dE/dx) from the MDC to calculate particle 
identification (PID) confidence levels (C.L.) for the z, K, 
and p hypotheses. The mass hypothesis with the highest 
PID C.L. is assigned to each track. A further requirement 
V, < 0.2 cm is imposed on the proton candidates to avoid 
backgrounds from beam interactions with residual gas in- 
side the beam pipe and materials of beam pipe and MDC 
inner wall. Photon candidates are reconstructed by clus- 
tering energy deposits in the EMC crystals. Good photon 
candidates are required to have energies larger than 25 
MeV in the barrel region (| cos 0| < 0.8) or 50 MeV in the 
endcap region (0.86 < | cos 0| < 0.92). To eliminate show- 
ers produced by charged particles, showers are required 
to be separated by more than 20? from anti-protons, and 
by more than 8? from other charged particles. The EMC 
time is required to be within (0, 700) ns of the event start 
time to suppress electronic noise and showers unrelated 
to the event [11]. 'The EMC shower shape variables are 
used to distinguish photons from anti-neutrons: the pho- 
ton candidates are required to have a lateral moment 
less than 0.4, and E3s3/ E55 larger than 0.85, where the 


Esx3 (Esx5) is the shower energies summed over 3 x 3 
(5 x 5) crystals around the center of the shower. 


In the studies of At > pr and At — pn? decays, the 
n mesons are reconstructed in their two most prominent 
decay modes, n > yy (Nyy) and n — m * m m? (n4 so), 
while the 7° meson is reconstructed in its dominant decay 
mode 7° — yy. Candidate  — yy and 4? — yy decays 
are selected using all yy combinations with an invariant 
mass within 3 times the mass resolution (10 (6) MeV/c? 
for the + (7°) signal) of their nominal masses (M, or 
Mo) id. An additional requirement, | cos Qaecay| < 0.9, 
where ñQaecay is the polar angle of one y in the helicity 
frame of the yy system, is imposed on the candidate 
7] — yy decay to suppress combinatorial backgrounds. 
To improve the momentum resolution, the yy invariant 
mass is then constrained to M, or M,o mass, and the 
resultant momenta are used in the subsequent analysis. 
The candidate n — mtmr’ are reconstructed using all 
TT T combinations with an invariant mass satisfying 
|Mr+r-70 — My| < 12MeV/c?. 

The Af is reconstructed using all combinations of the 
selected proton and the 7(7°) candidates. For et e~ anni- 
hilation at \/s = 4.6 GeV, there are no additional hadrons 
produced with the AZA- pair due to the limited phase 
space. Thus, two kinematic variables, the beam energy 
EP ag m FS jenes and 
the energy difference AE = Ext — Ebeam, are used to 
identify At candidates. Here, p at and E+ are the re- 
constructed momentum and energy of the AT candidate 
in the e*e^ c.m. system, and Eyeam is the energy of the 
electron and positron beams. For a A? candidate that is 
reconstructed correctly, Mpc and AE are expected to be 
consistent with the A? nominal mass and zero, respec- 
tively. A At candidate is accepted if the corresponding 
[AE] is less than 2.5 times its resolution (cA). The de- 
cay mode dependent AF requirements, are summarized 
in Table[I] For a given decay mode, we accept at most one 
charmed baryon candidate per event, retaining the one 
with the minimum |A E|. If there are candidates from 
different decay modes, we keep them all. For the decay 
mode At — pri44- 4», the peaking background from the 
CF decay mode A? > ntn X*(X* — pm?) is elimi- 
nated by requiring the invariant mass of proton and 7° 
satisfying |Mpro — My+| > 0.015GeV/c?. The MC study 
shows that the residual peaking backgrounds from At — 
ntr Xt (Xt — pr?) and from At — Artn’ (A — pr) 
and At > pK$n°(K3 — ntr), which have exactly the 
same final states as the signal, are negligible. 

The resultant Mpc distributions for the decays At — 
pn and At — pm? are depicted in Fig. [I] and Fig. 
respectively. The At — pr, signals are seen in both y de- 
cay modes, but no obvious A+ — pn? signal is observed. 
The data in the AE sideband region, defined as 3.50, g 
< |AE| < 6cAg, are used to study the backgrounds. 
The corresponding Mpc distributions, illustrated by the 
long-dashed histograms in Fig. [I]and Fig. P] show no At 
signal and that the combinatorial backgrounds are well 


constrained mass Mpgc = 
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described by the data in the AF sideband region. For 
the decay mode A? — pn,+,-7,0, data in the 7 sideband 
region (0.016 < |M,+,-70 — M,| < 0.032 GeV/c?), illus- 
trated by the (pink) dashed histogram in Fig. [I{b), also 
shows no evidence for peaking background. This is fur- 
ther validated by an analysis of the inclusive MC samples, 
where it is found that the combinatorial backgrounds are 
dominated by the processes e*e^ — qq. 
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FIG. 1. (color online) Simultaneous fit to the Mpc distri- 
butions of AZ — pn reconstructed with the decay modes (a) 
n — yy and (b) n > zT m). The dots with error bars are 
data, the (black) solid curves are for the best fits, the (blue) 
dash-dotted curves are for the backgrounds, and the (red) 
dashed curves are for the signals. The (green) long-dashed 
histograms and (pink) dashed histogram (in (b) only) are the 
data in the AE and M. -,o sideband region. 


To extract the signal yield for the decay At — pn, 
we perform unbinned maximum likelihood fits to the 
Mpc distributions. The signal probability density func- 
tion (PDF) is constructed by the signal MC simulated 
shape convoluted with a Gaussian function. Since MC 
simulation may be imperfect for modeling of the de- 
tector resolution and beam-energy spread of data, the 
mean and width of Gaussian function are free parame- 
ters to account for the potential mass shift and resolu- 
tion difference between data and MC simulation. The 
mean (ju) and width (c) values of Gaussian function are 
u = (0.74 + 0.56) MeV/c? and o = (0.32 + 2.28) MeV/c? 
for At — pny y, while u = (—1.22 + 0.80) MeV/c? and 
c = (0.02 + 1.44) MeV/c? for At > pr«4- 4o, respec- 
tively. The background shape is modeled by an AR- 
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FIG. 2. (color online) Fit to the Mgc distribution for the de- 
cay At — pr”. The dots with error bars are data, the (black) 
solid curve is for the best fit, and the (blue) dashed curve is 
for the background. The (green) long-dashed histogram is the 
data in the AE sideband region. The insert shows the nor- 
malized likelihood distribution, which includes the systematic 
uncertainty, as a function of the expected signal yield. The 
(blue) dashed arrow indicates the upper limit on the signal 
yield at 9096 C.L. 


GUS function [26] with the fixed high-end cutoff Epeam. 
The reliability of the ARGUS function is validated with 
the data in the AF sideband region as well as the in- 
clusive MC samples in the signal region. In the decay 
At — Pnr+r-70, the peaking backgrounds from the CF 
decays have been found to be negligible by MC studies, 
and are not considered in the fit. The fits are performed 
for the two 7 decay modes separately. The corresponding 
BFs are calculated using 


Naig 


B AT => E ; 
( c pn) 2 : NATA Mes Binter 


(1) 
where Nsig is the signal yield determined from the Mge 
fit, NATA = (105.9 + 4.8(stat.) + 0.5(syst.)) x 10? is 
the total number of Af A; pairs in the data [LI], e is 
the detection efficiency estimated by the MC simula- 
tion, and Binter is the 7 or n? decay BF taken from the 
PDG . The factor of 2 in the denominator accounts 
for the charge conjugation of the A+. Table[I]summarizes 
the signal yields, the statistical significances, estimated 
by the changes in the likelihood values obtained with and 
without the A? signal included, the detection efficiencies, 
and the resulting BFs. The two BFs for A? — pn, corre- 
sponding to the two 7 decay modes, are consistent within 
statistical uncertainties. 

We also perform a simultaneous fit to the Mgc distri- 
butions for the two 7 decay modes, constrained to the 
same B(A? — pn) and taking into account the different 
detection efficiencies and decay BFs of 7. The projec- 
tions of the fit curves are illustrated in Fig. [I] In the fit, 
the likelihood values of the two individual 7 decay modes 
are calculated as a function of BF, and are smeared by 
considering the correlated and uncorrelated systematic 
uncertainties (discussed in detail below) between the two 
n decay modes according to Refs. [2 Bs). The overall 
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TABLE I. Summary of the AE signal regions, the signal 
yields, the statistical significances, the detection efficiencies, 
and the BFs (where the first uncertainties are statistical, and 
the second systematic) for the different A7 decay modes. 


nn 


PNyy PNr+r— ro pm 
AE (GeV) [—0.034, 0.030] [—0.027,0.018]  [—0.056, 0.029] 
Nsig 38 +11 14+5 < 27.9 
Significance 3.20 2.70 n 
e(96) 39.8 20.3 49.0 
B(x10-3) 1.15 0.33 + 0.10 1.45 + 0.52 + 0.15 « 0.27 


likelihood value in the fit is the product of those for the 
two 7 decay modes. The resultant BF is determined to 
be B(A? — pn) = (1.24--0.28(stat.) +0.10(syst.)) x 1073 
with a statistical significance of 4.20, where the signifi- 
cance is estimated by the difference of maximum likeli- 
hood values for simultaneous fits with and without signal. 


Since no significant A+ — pma signal is observed, an 
upper limit on the BF is estimated. We fit the Mpa 
distribution for the candidate A? — pm? events using 
similar signal and background shapes to those described 
previously. The result of the best fit is shown in Fig. [2] 
For the signal PDF, the MC shape is convoluted with 
a Gaussian function with parameters fixed to those ob- 
tained in the fit to A? — pr, candidates. The PDF for 
the expected signal yield is taken to be the normalized 
likelihood £ obtained by scanning over the signal yield 
fixed from zero to a large number, and incorporating sys- 
tematic uncertainties Pg, as shown in the inset plot 
of Fig.[2] The upper limit at the 90% C.L. on the signal 
yield is NYP = 27.9 (shown as the arrow in Fig.[2), corre- 
sponding to m i £(x)dz/ J” £(z)dz = 0.9. The upper 
limit at the 90% C.L. on the BF is calculated with Eq. 
by substituting with 7° and is reported in Table [I] 

Several sources of systematic uncertainties are consid- 
ered in the BF measurements. The uncertainties as- 
sociated with the efficiencies of the tracking and PID 
for charged tracks are investigated with the samples 
ete” > (ntr), K* Kata and pprtn~ from data 
taken at \/s > 4.0 GeV, and the corresponding (trans- 
verse) momentum weighted values are assigned as the 
uncertainties. The uncertainties due to the V, require- 
ment and the veto on the CF peaking background in the 
decay At + pr44-4o are investigated by repeating the 
analysis with alternative requirements (V, < 0.25 cm and 
[Myzo — My+| > 0.020 GeV/c?). The resultant differ- 
ences of the BFs are taken as the systematic uncertain- 
ties. The z? reconstruction efficiency, including the pho- 
ton detection efficiency, is studied using a control sam- 
ple of D? + K-7*«? events from a data sample taken 
at ys = 3.773 GeV. The momentum weighted data-MC 
differences of the 7° reconstruction efficiencies, which are 
obtained to be 3.396 and 0.896 for At > pr44-4o and 
At — pm? decays, are considered as the uncertainties. 
Similarly, the uncertainty for the 7. reconstruction effi- 
ciency in the decay Af — pny. is determined to be 1.0% 


by assuming the same momentum-dependent data-MC 
differences as those for z? candidates. The uncertainties 
associated with the 7 mass window for At — pr 4— «o, 
the cos decay requirement for AT — pr44, the AE re- 
quirements, and the photon shower requirements are 
studied using double-tag D+ — m+m(z?) events. The 
uncertainties from the Mgc fit for At — pn candidates 
are studied by alternative fits with different signal shapes, 
background parameters, and fit ranges, and the resultant 
changes on the BFs are taken as the uncertainties. In the 
determination of the upper limit on the BF of At — pz? 
decay, similar alternative fits are investigated, and the 
one corresponding to the largest upper limit is selected 
conservatively. The uncertainties in the signal MC model 
arising from the following sources are considered: a) the 
beam energy spread; b) the input cross section line-shape 
of e*e^ — Mt A; production for ISR; c) the A? polar 
angle distribution in the efe rest frame; d) the differ- 
ent angular momentum between proton and 7(7°) candi- 
dates. The quadratic sum of the resultant differences in 
the detection efficiencies is taken as the uncertainty. The 
uncertainties of the MC statistics, the total AT Az num- 
ber quoted from Ref. and the decay BFs for the inter- 
mediate state decays quoted from the PDG are also 
considered. The total systematic uncertainties, quadratic 
sum of the individual ones, are 8.396, 10.296, and 5.296 
for AT > pnyy, priz-4-4o and pz), respectively. The in- 
dividual systematic uncertainties are summarized in the 
Table [I] 


TABLE II. Summary of the relative systematic uncertainties 
in percent for AZ > pnyy, pnz,«4-4,o and pr°. The sources 
tagged with ^*' are 100% correlated between the two 7 decay 
modes. 


Sources Dry  PNr+r-r0 pr 
* Tracking for p 1.3 1.3 1.3 
*PID for p 0.3 0.3 0.3 
Tracking for z*-— sine 2.0 ead 
PID for m+a7 w 2.0 T 
* V. requirement 0.2 0.2 0.2 
CF peaking background veto Se 1.3 ses 
nyy /n? reconstruction 1.0 3.3 0.8 
M +, zo mass window ET 1.2 

COS Ügecay requirement 1.2 - ias 
A E requirement 0.4 1.5 0.4 
Shower requirement 0.8 1.9 1.6 
Mpc fit 6.5 7.1 x 
Signal MC model 0.7 1.2 0.8 
MC statistics 0.1 0.1 0.1 
ai ATA 4.6 4.6 4.6 
Binter 0.5 1:2 negligible 
Total 8.3 10.2 5.2 


In summary, using 567 pb^! of e*e^ annihilation data 
taken at a c.m. energy of ys = 4.6 GeV with the BE- 
SIII detector, we find the first evidence for the SCS 
decay At — pn with a statistical significance of 4.20 
and measure its absolute BF to be B(A} — pn) = 
(1.24+0.28(stat.)+0.10(syst.)) x 1073. In a search for the 
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SCS decay At — pn, no obvious signal is observed and 
an upper limit at the 90% C.L. on its BF is determined 
to be B(A} — pz?) < 2.7 x 1074. The corresponding 
ratio of BFs between the two decays is also calculated 
to be B(A} — pr°)/B(At — pm) < 0.24, where the 
common uncertainties are cancelled. The measured BFs 
and their ratio are compared to the theoretical predic- 
tions from different models, as shown in Table [HI] Our 
measured BF of At — py is consistent, within two stan- 
dard deviations, with one of predictions in Ref. |3], the 
one that assumes flavor SU(3) symmetry and negative 
sign for p-wave amplitude of At > E°K*. It is worth 
noting that our measurement is significantly higher than 
other's theoretical predictions. The measured upper limit 
of B(A? — pn?) is compatible with the predicted values 
of most of theoretical models, but is smaller by a factor 
of 2 than that in Ref. Ep Overall, the obtained rela- 
tively large value of B(A% — pn) and the trend toward 
small value of the ratio B(A+ — pr°)/B(At — pn) will 
have a significant impact on theoretical calculation and 
will be helpful to understand the underlying dynamics 
of charmed baryon decays and to test SU(3) flavor sym- 
metry. Additional experimental data will improve the 
sensitivity of the measurements and allow a better dis- 
crimination between the different models. 


TABLE III. Comparison of measured BFs (in 107°) of AT > 
pn and pz? and their ratio to theoretical predictions. 


At pn Afp? quce 
Ad pn 
BESIII 1.24 + 0.29 « 0.27 « 0.24 
Sharma et al [3] 0.2*(1.7) 0.2 1.0*(0.1P) 
Uppal et al [4] 0.3 0.1-0.2 0.3-0.7 
S. L. Chen et al [12] e. 0.11-0.36° m 
Cai-Dian Lü et al [13] ks; 0.45 


* assume positive sign of p-wave amplitude of At > Z0 K+ 
assume negative sign of p-wave amplitude of At =O K+ 
© calculated relying on different values of parameters b and a 
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